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We consider synchronization of chaotic systems coupled indirectly through common environment 
where the environment has an intrinsic dynamics of its own modulated via feedback from the systems. 
We find that a rich variety of synchronization behavior, such as in-phase, antiphase, complete and 
anti- synchronization is possible. We present an approximate stability analysis for the difi'erent 
synchronization behaviors. The transitions to different states of synchronous behavior are analyzed 
in the parameter plane of coupling strengths by numerical studies for specific cases such as Rossler 
and Lorenz systems and are characterized using various indices such as correlation, average phase 
difference and Lyapunov exponents. The threshold condition obtained from numerical analysis is 
found to agree with that from the stability analysis. 

PACS numbers; 05.45.Xt 



U 



> 
00 

m 



On 
O 



X 



1. INTRODUCTION 

Chaotic synchronization of coupled nonlinear systems 
has been an area of intense research activity 1] . In such 
cases, depending on the strength and nature of coupling, 
the systems are capable of entering into different states 
of synchronization such as in-phase 0, H, antiphase [1| , 
lag [1, [Bl , anticipatory f/j, generalized sl-flO}. complete 
[ill ll2| and antisynchronization |13l-[l5 1 . Although all 
these different synchronization phenomena have been ex- 
plored in biological systems also, the case of phase syn- 
chronization is more useful in explaining many complex 
dynamical behavior in them. Specifically, antiphase syn- 
chronization with repulsive coupling has special relevance 
in biological systems such as neurons and ecological webs 

[IMl- 

Most of the present studies on synchronization con- 
sider mutually or unidirectionally coupled systems with 
or without parameter mismatch. However, synchroniza- 
tion has also been achieved by a common stochastic drive 
in uncoupled chaotic systems 0, 1131 • In such cases, 
the critical strength of noise for synchronization is nearly 
equal to the mean size of the attractor [2l| . The synchro- 
nized state thus often differs very much from the intrinsic 
characteristics of the individual system. Synchronization 
of chaotic systems by external periodic forcing where the 
driven system locks to the frequency of the drive has also 
been reported [l^-HJ- So also, a weak periodic force is 
found to stabilize inphase synchronization in a globally 
coupled array of Josephson i unctions [25j. 

Further, in the context of many real world sys- 
tems, synchronous behavior can occur due to interac- 



tion through a common medium. For instance, synchro- 
nization of chemical oscillations of catalyst-loaded reac- 
tants in a medium of catalyst free solution is reported 
where coupling is through exchange of chemicals with the 
surrounding medium [26| . So also, synchronized oscilla- 
tions in genetic oscillators occur due to coupling by diffu- 
sion of chemicals between cells and extracellular medium 
[27I [28j . Global oscillations of concentration of neuro- 
transmitter released by each cell can stimulate collec- 
tive rhythms in a population of circardian oscillators [2^ . 
Moreover, in an ensemble of cold atoms interacting with 
a coherent electromagnetic field, by controlling field cav- 
ity detuning, synchronized behavior with self-pulsating 
periodic and chaotic oscillations are found to occur [3(|. 
In all these cases, the coupling function has a dynamics 
modulated by the system dynamics. 

In general, such cases occur due to the common 
medium interacting with the dynamical systems. One 
refers to such a scheme as a coupling via a common en- 
vironment. The dynamics of n systems Xi, i = l,...,n 
coupled through an environment y is then given by 
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g{y) + h{xi,X2,.:,Xr. 
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where Xi and y have dimensions nix and ruy respectively. 
Such an indirect coupling has been reported in the con- 
text of periodic oscillators by Katriel [3l| . Under suitable 
conditions the periodic oscillators can synchronize. 

In this paper, we consider two chaotic systems coupled 
through a common dynamic environment as in Eq. ([l}. 
We show that this coupling can lead to a rich variety of 
synchronous behavior such as antiphase, in-phase, iden- 
tical, antisynchronization etc. This mechanism has the 
interesting feature that the common environment while 
capable of synchronizing the systems, does not cause ma- 
jor changes in their dynamics. In the synchronized state, 
the systems retain more or less the same phase space 
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structure of the uncoupled system. We present an ap- 
proximate stability analysis for the stability of the differ- 
ent synchronized states. We report detailed exploratory 
numerical studies for two standard systems, Rossler and 
Lorcnz, and demonstrate the rich synchronization behav- 
ior. The transition to different stages of synchronization 
is studied by computing average phase differences, cor- 
relations, and Lyapunov exponents. From the numerical 
studies, we verify the relation between the critical pa- 
rameters for the transition to different synchronization 
states obtained from the stability analysis. 



2. ENVIRONMENTAL COUPLING 

We consider two chaotic systems coupled to a common 
environment through a linear coupling 

ii = f{xi) + eilPiy (2a) 
±2 = f{x2) + eilP2y (2b) 
ij = -Ky- y7^(^ixi +/32a;2) (2c) 

The intrinsic dynamics of the environment is decaying 
with K as the damping parameter and therefore, without 
feedback from the systems, it is incapable of sustaining it- 
self for extended periods of time. Here, ei is the strength 
of feedback to the system and £2 that to the environ- 
ment. For simplicity, we take y to be one dimensional 
environment. Then, 7 is a column matrix {nix x 1), with 
elements zero or one, and it decides the components of 
Xi that take part in the coupling. 

The nature of feed back from and to the environment 
is adjusted by prescribing values for /3i and /32- When 
both fii and (32 are of the same sign, i.e. {Pi, (32) — 
(1, 1), the coupling is repulsive and can drive the systems 
to antiphase synchronization. When (3i and (^2 are of 
different signs, i.e. {(^1,(^2) = (1, —1), the coupling is of 
difference type leading to in-phase synchronization. We 
illustrate this behavior for the case of two chaotic Rossler 
systems coupled through environment as given by the 
equations 

in = -Xi2 - Xiz + ei(iiy 

ii2 = Xii + axi2 

ii?. b + x,:i{x,i -c) 

y = -K,y - y X! 

1=1,2 

The time series of the coupled Rossler systems for the 
in-phase synchronized and antiphase synchronized cases 
is shown in Fig. [T^ and[T]3. 

In the same way, two Lorenz systems are coupled 
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FIG. 1: Time series of the first variable xa of two envi- 
ronmentally coupled chaotic Rossler systems showing syn- 
chronization phenomena (a) In-phase synchronization (ei = 
£2 = 0.2, /3i = — /32 = 1) (b) antiphase synchronizatin 
(ei = £2 = 0.2, /3i = ^2 = 1). In both (a) and (b), we 
consider coupling only though one variable of the system, 
that is 71 = 1, 7i =0 for i ^ 1. Rossler parameters are 
a = b = 0.1, c = 18, i.e. we have chaotic attractor and the 
damping parameter, k = 1. 
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FIG. 2: Time series of the first variable xn of two en- 
vironmentally coupled chaotic Lorenz systems showing syn- 
chronization phenomena (a) In-phase synchronization (ei = 
£2 ~ 9.0, /3i = — /32 — 1) (b) antiphase synchronization 
(ei = 62 — 8.0, /3i = /32 = 1). Here, Lorenz parameters 
are (cr = 10, r = 28, 6 = 8/3). 



through environment as 

111 = a{xi2 - Xii) + ei(3^y 

112 = {r- Xi3)xu - xa 
xa ~ Xi\Xi2 — bxis 

y = -Ky - y X! ^^^'1 

1=1,2 

The in-phase and antiphase synchronized states of the 
coupled Lorenz systems are shown in Figs. [5^ and[2}D. 

This type of coupling is very relevant in the case of 
biological systems such as neurons where they interact 
through chemicals in the surrounding medium. We con- 
sider the case of two Hindmarsh-Rose model of neurons 
coupled through a common medium given by the follow- 
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FIG. 3: Time series of the first variable xn of two environ- 
mentally coupled Hindmarsh Rose neurons showing synchro- 
nization of bursts (a) In-phase synchronization (ei = £2 = 
0.4, l3i = — /32 ~ 1) (b) antiphase synchronization (ei — €2 — 
0.4, /3i = /32 = !)• Here, the parameters of the individual 
neuron are a = 3,6 = 5,r — 0.005,s — 4,c — 1.6,1 — 3.05 such 
that the individual neurons are chaotic. The synchronized 
state obtained is periodic in this case. 



ing equations 

il.2 = yi,2 + 0.x\ 2 - 2 - Zi,2 + I + eiPl,2W 



yi,2 = 1 - bxl ^ - yi,2 

zi,2 = -rzi,2 + sr {xi,2+ c) 

£2 \ ^ 



w 



-KW 



2 ^ 

1=1,2 



(^iXi 



(5) 



The in-phase and antiphase synchronized states of 
bursts for coupled HR neurons are shown in Fig. [5^ and 



3. LINEAR STABILITY ANALYSIS 

We analyze the stability of the synchronized state of 
two systems coupled via the scheme of Eq. If ^1, 

^2, and z represent the deviations from the synchronized 
state, their dynamics is governed by the linearized equa- 
tions obtained from Eqs. That is 

ii - /'(a;iKi+ei7/3i^ (6a) 

6 = I'{x2)^2+illp2Z (6b) 

z - -'t^- Y7^(/3ia+/326) (6c) 

In general, it is difficult to analyze the stability of the 
synchronized state from Eqs. ([6|). For the special case of 
the perfectly synchronized state, i.e. xi = X2, Eqs. ([H]) 
can be simplified by defining 

Co = Ml + M2- (7) 
Then Eqs. ([6]) can be written as 

e'o = f'{xi)Co+ei{pf+Pihz (8a) 
z = -nz- |7^eo (8b) 



The synchronized state corresponding to the fixed point 
(0, 0) of Eqs. ([8]) will be stable if all the Lyapunov expo- 
nents obtained from Eqs. ([8]) are negative. 

Considerable progress can be made if we assume that 
the time average values of f'{xi) and f'{x2) are approx- 
imately the same and can be replaced by an effective 
constant value A. In this approximation we treat and 
^2 as scalars. This type of approximation was used in 
Ref. [Ill and it was noted that it describes the overall 
features of the phase diagram reasonably well. Thus, us- 
ing Co defined by Eq. ([7]), Eqs. ^ can be written as 



Co = ACo + 2eiz 

(2^ 

where we choose j5\ + 0\ = 2. Eliminating z 
Eqs. (Pa]) and (|9b)) . we get an equation for Co as 

Co = (A - k)Co + (kA - eie2)Co 

Assuming a solution of the form 

Co = ^e"* 



(9a) 
(9b) 

from 
(10) 



we get 



(A - «) ± V(A - «;)2 - 4(6162 - A«;) 



(11) 



The synchronized state, defined by Co = /?iCi + (^2^,2 — 

0. is stable if Re[m] is negative for both the solutions. 
This gives the following criteria for the stability of the 
synchronized state. 

1. If (A — k)^ < 4(6162 — Ak), m is complex and the 
condition of stability is k > A. 

2. If (A — k)^ > 4(6162 — Ak), m is real and the stability 
condition becomes 6162 > Ak and k > A. 

In the first case above, the synchronized state is pos- 
sible if we have an environment which has a sufficiently 
fast decay to compensate for the divergence of the sys- 
tem due to A. In the second case, an additional condition 
must be satisfied. Here, the transition to stable synchro- 
nization is given by the threshold values of parameters 
satisfying the condition 



f2c 



Xk 



(12) 



We now consider the properties of the synchronized 
state defined by Co — /^iCi + /32C2 — 0, i.e. (3iXi + (^2X2 — 
const.. Numerical simulations show that the constant is 
zero. Thus, for /3i = /32 = 1 we get xi — —X2, i.e. an 
antiphase synchronization while for /3i = — /?2 = 1 we get 
xi = X2, i.e. an in-phase synchronization. 



4. NUMERICAL ANALYSIS 

The scheme of coupling through the environment given 
in Eqs. © is applied to standard Rossler and Lorenz 



systems. We study the two cases, /3i = +1 and /32 = —1 
where in-phase synchronization is possible and Pi — 1^2 = 
+1 where antiphase synchronization is possible. 

When /3i = +1 and /32 = —1, we observe in-phase 
synchronization in both Rossler and Lorenz systems 
(Figs, [ija) and[2][a)). As the coupling strength is in- 
creased, systems go to a state of complete synchroniza- 
tion. When /?! = /32 = +1 the synchronized states are 
out of phase with each other giving antiphase synchro- 
nization for both Rossler and Lorenz (Figs. [IJb) and 
^b)). As the strength of feedback is increased in the 
case of Rossler systems, control of chaos is observed and 
the systems become periodic, but the two coupled sys- 
tems are still in antiphase synchronization. In the case 
of Lorenz systems as the coupling strength is increased, 
the systems become anti-synchronized where xi — —X2, 
yi = -2/2 and 2:1 = 22. 



4.1. In-phase (or Antiphase) synchronization 
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FIG. 4: Transition from regions of no synchronization to 
antiphase synchronization is shown in the parameter plane 
ei — £2 for the coupled Rossler systems. The points are ob- 
tained numerically when the phase difference becomes approx- 
imately TT. Solid curve corresponds to the stability condition 
Eq. dill, i.e. e2c oc 1/eic. 



The transitions to in-phase ( or antiphase ) synchro- 
nization can be studied numerically using the average 
phase difference between the two systems. For this, we 
need to define phases of individual systems. In the case of 
Rossler systems, as the trajectory has a rotation around 
a fixed point in the x — y plane, the phase ^(i) of the 
Rossler system can be defined [3J| as the angle 



0(i) =tan-i(j/(t)/x(i)) 



(13) 



The phase (pit) and the phase difference ipit) between 
the two Rossler systems coupled through environment 
are calculated using Eq. (IT3l) for increasing strengths of 
feedback for identical feedback strengths ei ~ €2- The 
mean phase difference over many cycles < V'(^) >j is 
nearly for the in-phase synchronization and tt for the 
antiphase synchronization. 

Since the Lorenz system does not have such a proper 
rotation around any fixed point, the phase cannot be 
defined by Eq. ((T3)) . The phase of Lorenz system is cal- 
culated using the modified variables [2^ as 



(f){t) — tan ^{z/u) 



(14) 



where u = u- 
p—1 and u = 



V2/3(p - 1), 2p = 
^ The dynamics in (w, z) looks like 



and Up 



a rotation around some center point (up, Zp). The phase 
(j){t) of the individual Lorenz systems are calculated using 
Eq. (fT4|) . The phases show confinement due to coupling 
indicating in-phase ( or antiphase ) synchronization. It 
is evident that since we neglect the sign of x and y in 
the calculation of m, phase defined as in Eq. can 
not distinguish between in-phase and antiphase cases. In 
this context, the similarity function S [H and a modified 
similarity function S' [35| serves as a useful index for 
identifying the in-phase or antiphase synchronization. 



The similarity function S is defined for a delay time r 

, _ < [X2{t + t) - Xi{t)Y > . . 

- [< Xl{t)><xl{t) >]V2 (15) 

and the modified similarity function S' is defined as 



S\r) 



< [X2{t + T) +Xl{t)]^ > 
[< Xlit)><xlit) >]V2 



(16) 



For Pi = 1,P2 = —1, at T = 0, S = corresponding 
to the complete synchronization and S is finite for the 
in-phase synchronization. Similarly, for Pi — P2 = 1, &t 
r = 0, S" is indicating the antisynchronization and S" is 
finite for the antiphase synchronization. For the coupled 
Rossler systems the average phase difference is calculated 
for the full parameter plane (£1,62) in the range (0,0.5) 
and the points where the value becomes approximately 
TT is plotted in Fig. 2] These therefore correspond to 
the threshold values for onset of stability of antiphase 
synchronizaton. The full line corresponds to the curve 
plotted using the threshold condition from our stability 
theory in Eq. (1121) . The agreement is quite good with a 
A = 0.009 and the relation e2c oc 1/eic is clearly seen. 
Similar transition curves are also observed for transition 
to in-phase synchronization in the case of /3i = +1 and 
P2 — —1 and also for Lorenz systems and they agree 
with the relation e2c oc 1/eic obtained from the stability 
theory. 

As seen in Eq. (fT2| . we also have the relations €2c oc k 
and eic oc k. Figure. [5] shows the phase plot for the transi- 
tion from unsynchronized to antiphase synchronized state 
in the ei — k plane. A linear relation is clearly seen and 
the solid line is drawn with the effective A = 0.009, thus 
validating the transition criterion of Eq. (|12l) obtained 
from the stability theory. 
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FIG. 5: Transition from unsynchronized to antiphase syn- 
chronized regions is shown in the parameter plane (k, ei) for 
coupled Rossler systems. Points are obtained from numeri- 
cal simulation with £2 = 0.6 and the solid curve is a linear 
fit corresponding to the stability condition Eq. (|12p with the 
effective A = 0.009. 



4.2. Lyapunov exponents 
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FIG. 6: Four largest Lyapunov exponents are shown as a 
function of the feed back strength e for the two coupled 
Rossler systems coupled through dynamic environment with 
ei — €2 = e. (a) Pi — +I,l32 = —1; the first crossing of 
at e = 0.12 (fourth largest LE) indicates the transition 
to in-phase synchronization, while the second zero crossing 
at e = 0.59 (second largest LE) indicates the transition to 
complete synchronization (b) /3i = P2 = +1; the first cross- 
ing of at e = 0.12 indicates antiphase synchronization and 
the region where all Lyapunov exponents are less than or 
equal to zero indicates the antiphase synchronized periodic 
states. (Lyapunov exponents are calculated by considering 
variational equations using Wolf algorithm [s^l-) 



The transitions to all the different types of synchro- 
nization described above can be tracked by calculating 
the Lyapunov exponents. Since the coupling here is in- 
direct and through an environment, instead of calculat- 
ing transverse Lyapunov exponents about the synchro- 
nized state, we calculate all the Lyapunov exponents of 
the coupled system. The variation of these Lyapunov 
exponents with coupling strength helps to identify the 
onset of in-phase ( or antiphase ) and complete ( or 
anti) synchronization. The two chaotic systems and the 
environment together form a seven-dimensional system. 
The changes in the largest four Lyapunov exponents are 
used to identify transitions to different synchronization 
regimes. First crossing from zero to negative of the fourth 
Lyapunov exponent indicates the onset of in-phase ( or 
antiphase ) synchronization and the crossing of the sec- 
ond largest Lyapunov exponent indicates the onset of 
complete ( or anti ) synchronization The largest 
four Lyapunov exponents for coupled Rossler systems are 
shown in Fig. |6] for various strengths of feedback. For 
the case (3i — +1 and /32 — —1, the zero crossing of 
the fourth largest Lyapunov exponent in Fig. El^a) corre- 
sponds to the onset of in-phase synchronization, and the 
zero crossing of the second largest Lyapunov exponent 
corresponds to the onset of complete synchronization. 
Here, the narrow window where all Lyapunov exponents 
are less than or equal to zero corresponds to synchro- 
nized periodic states in Rossler systems as verified from 
the time series. In the case of antiphase synchronization 
similar results are seen (Fig. IH^b)). The region where 
all Lyapunov expoents are less than or equal to zero in 
Fig. ^h) corresponds to the antiphase synchronization 
in the periodic state. The results of a similar analysis 
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FIG. 7: Four largest Lyapunov exponents are shown as a 
function of the feedback strengths e of two coupled Lorenz 
systems with ei = £2 = e. (a) /3i = +I,p2 = —1; the first zero 
crossing of Lyapunov exponent at e = 4.2 indicates in-phase 
synchronization, the second zero crossing at e = 9.8 indicates 
complete synchronization (b) /3i = (32 — +1; the first zero 
crossing at e = 4.4 indicates antiphase synchronization, the 
second zero crossing at e = 9.8 indicates antisynchronization. 

for Lorenz are shown in Fig. [T] In Fig. El^a), the case 
/?i = +1 and = — 1 is shown where the zero crossing 
of the fourth largest Lyapunov exponent corresponds to 
the onset of in-phase synchronization, and the zero cross- 
ing of the second largest Lyapunov exponent corresponds 
to the onset of complete synchronization. In Fig. [Tljb), 
the case /3i = /32 = 1 is shown where the zero crossing of 
the fourth largest Lyapunov exponent corresponds to the 
onset of antiphase synchronization, and the zero crossing 
of the second largest Lyapunov exponent corresponds to 
the onset of antisynchronization. 




4.3. Phase diagram 

In this section, we present the complete phase diagram 
in the parameter plane of coupling strengths identifying 
the regions of different states of synchronization such as 
complete ( or anti ) synchronization, in-phase ( or an- 
tiphase ) synchronization and unsynchronized regions. 
We use the average phase difference and Lyapunov ex- 
ponent to mark the different regions of synchronization. 
In addition, the complete and antisynchronization states 
are characterized by calculating correltion between the 
two systems using 

^ ^ (Mt)- < Xi{t) >){x2{t)- < X2{t) >) 

The phase diagram in the ei - £2 plane for Rossler sys- 
tem is shown in Fig. for /3i = 1, /32 = — 1- As 
the coupling strengths increase ( along the diagonal ) 
we see a transition from the unsynchronized state (dark 
gray) to the in-phase synchronized state (light gray) and 
then to the completely synchronized state (white). For 
large coupling constants, the system becomes unstable 
(black). The critical coupling constants corresponding 
to the transitions between the different types of synchro- 
nization obey the relation Eq. ([T^ as obtained from the 
stability analysis. Figure. [8J5 shows a similar phase dia- 
gram for pi = /32 — 1. Here dark gray region corresponds 
to unsynchronized states, region marked I corresponds 
to antiphase synchronization in chaotic state, regions II 
- IV corresponds to different regimes of synchronization 
in periodic states and black region corresponds to unsta- 
ble states. We find that here, depending on the coupling 
strength the coupled systems settle to two different peri- 
odic states A and B. The x — y plane corresponding to 
the states A and B are shown in Fig. [5] In regions II and 
IV, both are in state A shown in Fig. [3^. while in region 
III, one system is in periodic state A and the other in 
state B{Fig.Mi). 




FIG. 8: Regions of different states of synchronization marked 
out in the parameter plane (ei — €2) for the coupled Rossler 
systems. The different phase space regions are obtained by 
using the asymptotic correlation values, average phase differ- 
ences and Lyapunov exponents, (a) Pi = +1,(52 = —1. White 
region corresponds to |C| ~ 0.99 indicating synchronized re- 
gions; light gray region is in-phase synchronized region, (b) 
Pi = 132 = +1 region I corresponds to antiphase synchronized 
chaotic states, regions II - IV corresponds to different states 
of antiphase synchronized periodic states (see text). In both 
cases, the dark gray region corresponds to the unsynchronized 
states and the black regions in the upper right corner are the 
unstable states. 
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FIG. 9: The x — y phase plane of antiphase synchronized 
periodic states in regions II and III of Fig. |HJd. (a) ei = £2 = 
1.5 Both systems are in state A (b) ei = £2 = 2.5 Systems are 
in different states A and B. 



In regions II and IV, the synchronized states are such 
that xi{t -I- r) ~ X2{t), corresponding to lag synchro- 
nization and in region III, the systems are in antiphase 
synchronization in the periodic state. The average error 
function calculated after shifting xi{t) by half the time 
period for the regions I, III and IV is shown in Fig. [TUl 



The similar phase diagrams for coupled Lorenz sys- 
tems are shown in Figs. [TTb and [Tib . As e is increased 
along the diagonal, we observe transitions in the follow- 
ing sequence: unsynchronized state (dark gray) to in- 
phase/ antiphase synchronized states (light gray) to com- 
plete/anti synchronized states (white) to unstable states 



(black). Here also, the critical coupling constants corre- 
sponding to the transitions between the different types 
of synchronization obey the theoretical relation Eq. (|12l) . 



5. CONCLUSION 

We report the synchronization of two nonlinear chaotic 
systems by coupling them indirectly through a common 
environment. The coupling mechanism proposed is gen- 
eral and can be adjusted for in-phase and antiphase or 
complete and anti- types of synchronization. The dif- 
ferent types of synchronous behavior and the transitions 
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among them are analyzed in the case of two standard sys- 
tems Rossler and Lorenz using the numerically computed 
Lyapunov exponents, average phase difference, correla- 
tion from time series and similarity function. Using an 
approximate linear stability analysis, the threshold val- 
ues of coupling strengths for onset of synchronization of 
the in-phase or antiphase type are derived and the tran- 
sitions curves obtained from numerical calculations agree 
with the curves from stability analysis. 



FIG. 10: Average error function computed after shifting one 
of the time series by half the period in the synchronized pe- 
riodic regions II - IV of Fig. [Sjs. The average error ~ for 
e < 2 and e > 2.95 indicating lag synchronization. The region 
2 < e < 2.95 corresponds to antiphase synchronization in the 
periodic state. 
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This method of synchronization has the interesting fea- 
ture that the synchronized state has almost the same 
phase space structure as that of the uncoupled dynamics. 
The mehod reported here offers a simple coupling scheme 
to realize phase ( or antiphase ) synchronization in two 
coupled chaotic identical systems. As far as we know, the 
reported works in this are mostly on nonidentical systems 
with parameter mismatch or delay in coupling. 



FIG. 11; Regions of different states of synchronization marked 
out in parameter plane of coupling strengths (ei, €2) by com- 
puting asymptotic correlation values, average phase differ- 
ence and similarity functions for Lorenz systems, (a) jSi — 
+ I,j32 = —1 (b) Pi — l32 = +1. In both cases, white re- 
gion corresponds to |Cj ~ 0.99 indicating synchronized/anti- 
synchronized regions. Light gray region is in-phase ( or an- 
tiphase ) syncrhonized state and dark gray region is unsyn- 
chronized state. Black region in the upper right corner cor- 
responds to unstable states. In the numerical simulations, ei 
and £2 are varied in steps of 0.2. The time averages in Eq. (|17p 
are taken over 50 time units. 



The results for synchronized states with such a cou- 
pling are presented here for three standard cases such 
as Rossler, Lorenz and Hindmarsh-Rose systems. How- 
ever, we have checked that it works in general for a few 
cases also such as FitzHugh Nagumo model of neurons 
and Mackey- Glass system. 
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